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,^ Abstract: Despite being very successful in explaining the wide range of precision exper- 

imental results obtained so far, the Standard Model (SM) of elementary particles fails to 
^]^ address two of the greatest observations of the recent decades: tiny but nonzero neutrino 

(^ masses and the well-known problem of missing mass in the Universe. Typically the new 

^^ models beyond the SM explain only one of these observations. Instead, in the present arti- 

1/-^ cle, we take the view that they both point towards the same new extension of the Standard 

ly-j Model. The new particles introduced are responsible simultaneously for neutrino masses 

^^ and for the dark matter of the Universe. The stability of dark matter and the smallness 

,__l of neutrino masses are guaranteed by a U(l) global symmetry, broken to a remnant Z2. 

^ The canonical seesaw mechanism is forbidden and neutrino masses emerge at the loop 

'^ level being further suppressed by the small explicit breaking of the U(l) symmetry. The 

new particles and interactions are invoked at the electroweak scale and lead to rich phe- 
nomenology in colliders, in lepton flavour violating rare decays and in direct and indirect 
dark matter searches, making the model testable in the coming future. 

Keywords: Beyond Standard Model, Neutrino Physics, Cosmology of Theories beyond 
the SM. 
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1. Introduction 

Despite all its triumphs, the Standard Model (SM) of elementary particles fails to explain 
two of the greatest observations of the recent decades: tiny but nonzero neutrino masses 
and the missing mass of the Universe commonly explained by Dark Matter (DM). Various 
models beyond the Standard Model (BSM) have been developed to explain each of these 
mysteries separately. On the contrary these two phenomena may be linked and explained 
within a single scenario [1, 2, 3]. In the present article we will explore this possibility and 
propose an extension of the Standard Model in which the new particles and interactions are 
simultaneously responsible for the dark matter of the Universe and neutrino masses. We 
call our model AMEND which stands for "A Model Explaining Neutrino masses and Dark 
matter". We introduce a \]{l)x global symmetry which is broken to a remnant Z2. The 
latter symmetry distinguishes the new particles from the SM ones and is responsible for 
the stability of dark matter. At the same time, it forbids the canonical seesaw mechanism 



and left-handed neutrinos do not acquire a Dirac mass. Neutrino masses arise at the loop 
level and are suppressed by small terms which explicitly break the U(l)x symmetry. 

More concretely, we introduce two additional electroweak doublets R and R' with 
opposite hypercharges. These Weyl fermions together can be regarded as a Dirac four 
component spinor. Moreover, we add an electroweak scalar triplet A and a complex singlet 
4> which do not acquire vacuum expectation values (VEVs). The lightest neutral scalar, 
which is mainly the singlet (p with a small admixture of the neutral component of A, 
will play the role of the DM particle. We expect the mass value for DM to be around the 
electroweak scale so the DM in our model can be categorised as a weakly interacting massive 
particle (WIMP). Smaller values can be obtained if fine-tuning is allowed in the DM mass 
term. The Yukawa coupling involving A is the lepton number violating interaction which 
leads to Majorana neutrino masses. The neutrino masses are generated at one loop and are 
proportional to the mass of the new fermionic doublet which needs to be at the electroweak 
scale or higher. The model presented in this paper belongs to the class of models for which 
the neutrino mass is generated radiatively (see e.g. [4]). The low neutrino mass scale 
therefore requires a further suppression, in addition to the loop-factor, by either suppressing 
the Yukawa coupling or by requiring a cancellation between the various contributions. In 
our model both of these requirements are simultaneously enforced by a continuous Abelian 
symmetry U(l)x which forbids the neutrino masses. Its explicit breaking to a residual Z2 
symmetry leads to small neutrino masses and a stable DM particle at the same time. 

In our model, DM is produced thermally in the Early Universe. We explore the various 
allowed annihilation channels which can proceed via R-R' , Z or Higgs exchange, in order 
to reproduce the observed amount of dark matter r^oM^^ = 0.1131 it 0.0034 [5]. We find 
that in the parameter range of interest, the dark matter annihilation through s-channel 
Higgs exchange dominates over the other channels. Dark matter is being searched for both 
directly, looking for the recoil of nucleons/electrons due to DM scattering in detectors, and 
indirectly, observing the products of DM annihilations (photons, positrons, neutrinos, anti- 
protons, anti-deuteron) in overdense regions in the galaxy or inside astrophysical objects 
such as the Sun or the Earth. No positive signal has been found in direct searches except 
for the two events in CDMS-II [6] and the DAMA experiment which has reported a positive 
signal at 8.2a [7]. Very recently, the CoGeNT collaboration [8] has found an indication 
for excess of events that might be due to scattering of light DM off the nuclei. The first 
analysis of the XENONIOO experiment [9] excludes all positive DM signals, however the 
interpretation of the data depends on the astrophysical uncertainties [10] as well as the 
effective light yield in the low mass region [11]. We will consider these possible signals 
for DM detection and check whether they can be accommodated in our model. In the 
case of DAMA, we look into two possible explanations, by either light dark matter or 
by the inelastic scattering scenario. As our candidate is a quasi-singlet coupled to the 
Higgs boson, the interactions which induce the dark matter annihilations at freeze-out are 
also responsible for elastic or inelastic scattering off nuclei, relevant for direct dark matter 
searches. In some cases the cross sections of the two processes, annihilations and scattering 
off nuclei, are related. 

All new particles are typically expected to have masses around the electroweak scale 
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Table 1: Particle content and gauge quantum numbers. 



and to couple to the SM particles at tree level leading to rich phenomenology which makes 
the model testable in the near future. We therefore investigate several experimental bounds 
including the invisible decay width of the Z into DM pair (if kinematically allowed) , the 
branching ratios of Lepton Flavour Violating (LFV) rare decays and the anomalous mag- 
netic moment of the muon. Furthermore, we comment on the possible signatures at the 
LHC. It is possible to correlate the flavour structure of the couplings measured at the LHC 
with the neutrino mass matrix and the data from LFV rare decay searches. 

The paper is organised as follows. In sec. 2, the model is presented. In sec. 3, the 
neutral scalar sector is analysed. In sec. 4, the neutrino mass generation at one loop level 
as well as effects on the LFV rare decays and magnetic dipole moment of the muon are 
discussed. In sec. 5, different processes that can give rise to annihilation of dark matter are 
explored. A discussion of the possibilities of direct and indirect dark matter detection is 
also included. In sec. 6, experimental constraints from electroweak precision tests as well 
as possible collider signatures are studied. Finally, in sec. 7, results are summarised. 

2. The Model 

In order to explain neutrino masses and dark matter, we extend the SM with two additional 
scalar fields, and one vector-like fermionic doublet. The complete particle content of the 
model and the SM quantum numbers are summarised in Tab. 1. More specifically, the 
scalar sector of the model contains three fields: 



the SM Higgs doublet which is indicated by H in the following; 

a complex field, (j) = {4>i + i02)/\/2j which is a singlet of 'S\]{2)l x U(l)y; 

and a triplet scalar field A: 



(2.1) 



where the neutral component can be decomposed as A'^ = (Ai + i A2)/v2, with Aj 
being real fields. 

In the fermionic sector, the added vector-like SU(2)i doublet is described by two Weyl 
fermion SU(2)l doublets, R^ = {vr E]^) and {R')'^ = {E+ z/^). 

With this particle content, a model enjoying a very high level of symmetry can be 
constructed. We consider a Lagrangian which preserves the SM gauge group as well as 
U(l)£ of lepton number, \]{1)^ under which only (j) is charged, a similar U(1)a for A and 
U(l)/j under which R and R! have opposite quantum numbers. Let us define 

G = U(1),jxU(1)<^xU(1)axU(1), . (2.2) 

The G-preserving part of the scalar potential is given by 

y = - fij^H^H + /iitr (At a) + fil^U 

+ ^{H^Hf + ^{4>Uf + — (trAtA)' + ^tr(At[At, A]A) 
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+ XHAiH^Hti (At a) + AHA2^"^[At , A]H + A<^a0V tr (AtA) + \h4>c^Hh^ H , 

(2.3) 

and the fermionic part contains the Dirac mass term of the vector-like doublet 

-^R = ruRRiR'^)^ • R + h.c. , (2.4) 

where {R ) = {v'^ — (E^) ). In order to avoid present collider bounds, we require 
1TT-RR to be larger than ~ 100 GeV. For definiteness we will take itirr^ = 300 GeV, unless 
otherwise stated. Terms in Eqs. (2.3) and (2.4) constitute the most general renormalisable 
gauge invariant Lagrangian preserving G that can be added to the SM Lagrangian. 

Among all possible U(l) subgroups of G which can be obtained by assigning different 
possible charges to the fields, we list a number of symmetries that are of particular interest 
in Tab. 2. Notice that the quarks and the SM Higgs field have zero quantum numbers 
under these symmetries. We assume a hierarchical pattern for the breaking of the group. 
First at a very high energy, Ah, the group G breaks to U(l)x under which the fields are 
charged as in Tab. 2. Note that Gsm x U(1)x is anomaly-free as the new fermionic doublet 
is vector-like. The terms which arise after G — )■ U(l)x are 

^HA^ =XHA<pH^icJ2A^H(t>^ + h.c. (2.5a) 

-^i^^P =ga<t>^R^iLa + h.c. . (2.5b) 

After electroweak symmetry breaking, the first term will induce mixing between (p and A. 
The second term introduces a coupling between the new sector and the leptonic doublet. 
This U(l)x symmetry is eventually broken into a residual Z2, under which SM particles 
are even and the new states are odd. The Z2, being exact, forbids a Dirac mass term of form 
R^iia for neutrinos and neutrino masses cannot therefore arise from the seesaw mechanism. 



Moreover, it guarantees the stability of the hghtest new particle which is a potential dark 
matter candidate. At low energy, we expect a theory which is nearly U(l)x-conserving 
with small breaking terms which preserve Z2. We assume this breaking to be explicit for 
the purpose of the present study. Notice that the terms ^^0 and YnAcf) respect U(l)x and 
are not therefore suppressed. The U(l)x-violating contributions to the scalar potential are 

■Scalar = ~XHA<pH^ia2^^ H(t> + fil<f>'^ + \^i<f>'^ + \^24>^4>^ + \h^H^ H^^ +~XA<ptTA^ A^^ + h.C. . 

(2.6) 
The new Weyl fermions R and R' couple to the SM leptons with two additional Z2- 
preserving terms 

-^e^^P = ga(pR^iLa + h.C. and - M^A = {gA)aR'^ ■A-iLa+ h.c. . (2.7) 



Due to the assumed breaking pattern of the G symmetry, we have the hierarchy g ^ g,gA 
and Xha4> ^ XhA(P- The freedom of a global phase transformation of and A can be used 
to set the phases of \HA<f> and /i| to zero. Moreover, the phases of ga can in general 
be absorbed by (.La- Thus, the U(l)x-preserving part as well as the mass terms can be 
made real. In this basis, cja and {gA)a can in general be complex leading to CP-violating 
Majorana and Dirac phases in the neutrino mass matrix. In this paper, for simplicity we 
restrict our analysis to the CP conserving case. Notice that, in general, the couplings X^i, 
X(f,2, Xh4>-, Xa4>, Xha4> can be either positive or negative. The couplings have to be taken in 
a range such that the potential is stable at infinity. Since the A^i, A^2) Xutp-, Aa0 and \ha<j> 
couplings are much smaller than the corresponding U(l)x-conserving terms, the potential 
remains stable regardless of their sign. Unless otherwise specified, we take these couplings 
to be positive in our studies. A similar analysis and similar results could be obtained for 
negative couplings. 



3. Neutral Scalar Masses 

The parameters of the model can be chosen such that only the SM Higgs field develops a 
vacuum expectation value, that is 
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Table 2: Specific U(l) sub-groups of G and associated particle quantum numbers. Indicated are 
also the terms in the full Lagrangian which violate each symmetry. 



As a result, the Z2 symmetry is preserved. The Z2 symmetry prevents mixing between 
the new scalars and the SM Higgs but the Xha^ and XhAcj) couphngs in yHA<j> and yH£:^(i> 
lead to mixing between (/> and the neutral component of A. There are four massive neutral 
scalar fields in the model, (5i_2,3,4) with masses respectively given by 
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For simplicity, we have assumed ml — tti^ ^ ^'SA ^^ ^^^^ ^^ "^A ^ "''i- Notice also 
that 771^ » m^ and w,^^ ^> m^^, as the parameters indicated by tilde are the U(l)x- 
breaking ones. Thus, the pair of states {5i, 82) and {63, 5^) are nearly degenerate with 
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We define the mass splitting 
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The mass eigenstates are 
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where, at leading order, | tan2ai| ~ | tan2Q:2| — 2m^^/(mA ~ '^I)- The difference in |ai| 
and \a2\ is suppressed by the U(l)x-breaking terms. Notice that, if the couplings Xhi/, and 
XhA(j) were taken to be negative, the roles of 5i and 82 as well as those of (Js and 84, would 
be interchanged, with 82 being the lightest particle and the dark matter candidate. In 
appendix A, we describe in detail the mass matrix and give the general mass eigenvalues of 
the scalars and the mixing between {(j)i, 02, Ai, A2) with the mass eigenstates {81,82-, (^s, 8^). 
By taking //^ relatively large, the components of the triplet can be sufficiently heavy and 
the bounds from direct searches can be therefore avoided. In our model, the values of Mi 



and M2 are considered free parameters and can range from a few keV (in order to avoid 
too hot dark matter) to above the electroweak breaking scale. 

For hght scalar masses with Mi + M2 < mz, there will be an additional invisible decay 
mode of the Z boson into 6162 due to the coupling of A to the Z boson 

i 9sTi('9'i sin ai sin qo 

^^ ^ [628^,5^ - 6,d^52]Z'^ , (3.6) 

where gsu{2) is the SM weak gauge coupling and 6w is the Weinberg angle. The corre- 
sponding decay width is given by 

Gf sin^ ai sin^ 02 3 . . 

T{Z -^ 6162) = -^= m| , (3.7) 

which is strongly sensitive to the mixing between (j) and A, i.e. sinai sin 02. 62 eventually 
decays into 5i and neutrinos via Z or R — R' exchange as it is given in Eq. (5.17). The whole 
process appears as a Z invisible decay mode. Hence the present bound on the invisible 
decay width [12] constrains the mixing sinoi sin 02 as 

T{Z -^ 8182) < 0.3% Tin^isibic => sinai sin 02 < 0.07 . (3.8) 

For heavier masses, this bound does not apply and larger mixing is in principle allowed. 

For definiteness in the following study we will take the following typical values for the 
scalar parameters, unless otherwise indicated: 

Ml ~ M2 = 70 GeV , M3 ~ M4 ~ mA = 500 GeV , <5 = 50 MeV , 

sinai ~ —sin 02 = —0.1 . (3.9) 

In principle, in our model, the DM particle can have much smaller masses if strong fine- 
tuning is allowed in the masses (see Eq. (A.4a)). In this case, the lower bound on the mass 
is given by large scale structure formation, i.e. few keV (see e.g. [13]), and by big bang 
nucleosynthesis, depending on the dominant DM couplings to SM particles [14] . We do not 
consider further this case in our study and we focus only on electroweak-scale DM masses. 

4. Lepton Sector 

4.1 Neutrino Masses 

Neutrino masses are protected by the symmetry of the model. If U(l)/j x U(l)(^ x U(1)a x 
11(1)^ is conserved, neutrinos cannot possess a Majorana mass term. The unbroken Z2 
symmetry prevents a Dirac mass term such as R'iLa- The term in ^e^^f, is allowed by 
U(l)x but does not generate a neutrino mass term as neither (p nor A acquire a vacuum 
expectation value, {(/)) = (A) = 0. We notice that if either of U(1)li, U(1)l2 or \J{1)l3 
(defined in Tab. 2) were conserved, the neutrino mass would be protected by a lepton num- 
ber symmetry. Once the symmetry is explicitly broken, a Majorana mass term can emerge. 
We therefore expect the neutrino mass to depend on combinations gA9^HA4ii dAd^HAcp^ 
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Figure 1: Effective neutrino mass generation at one loop. 



dAg^HA^ and gAg^HA(t>- Thus, the smahness of neutrino masses can be explained by 
t'Hooft's criterion [15]. 

More specificahy, an effective neutrino mass term, 

-^u^u^ = 2 ("^-)a/3 {^l)a ^ (^i)/3 + ^-C , (4.1) 

arises at one loop- level [4, 1, 2, 3] through the diagram shown in Fig. 1. The neutrino mass 
matrix is given by 

{m^)al3 = [ga{gA)/3 + gfi{gA)a\fi + [ga{gA)fS + gfi{gA)a\ri , (4.2) 

and depends on g^ as expected. The determinant of the neutrino mass matrix vanishes, 
so one of the mass eigenvalues is zero, unless more vector-like fermionic doublets or copies 
of (/) or A are added. In other words, within the present model with only one generation 
of R and R' fields, the neutrino mass spectrum is either normal hierarchical or inverted 
hierarchical. The terms r/ and fj can be explicitly computed 

mBR ( Mi m\^^ Ml m\^^ 
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It is straightforward to check that in the limit in which \ha4> and Xha<i> both vanish, i.e. 
m^^ = fTij,^ = 0, the neutrino mass becomes zero. This is expected as in this case the 
symmetry U(l)i3 is exact (see Tab. 1 for definition) and no Majorana mass for neutrinos is 
allowed. In the limit of nearly-exact U(l)x symmetry, we have Mf ~ M^ ~ "^i~"^lA/''^A' 



M^ ~ M^ and sinai ~ — sin 02- As a result, in this limit we can approximately write 
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We have expanded to first order in the U(l)x-breaking parameters and assumed that 



m 



i^ <^ mj^ — mi. In Eqs. (4.4b) and (4.4d), we have taken the limit mi, Mf <C mj^, ^^rr 
in addition. Notice that r) is suppressed by U(l)x-violating parameters 'fni^/ {m^^j^^ — m,\) 
and 'fhi/m?^^ relative to r/. On the other hand, the contribution of rj to the neutrino 
mass is suppressed by the U(l)x-violating coupling g <^ g. As expected, the neutrino 
mass depends on the mixing between cj) and A (given by the terms XhAcP and XHAtp after 
electroweak symmetry breaking) and on the coupling between the new sector with the 
leptonic doublet (i.e., g, g and ^a)- We can obtain an order of magnitude estimate for the 
couplings 
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where we have taken as typical values m.RR = 300 GeV and m.^, = 0.05 eV in Eqs. (4.4b) 
and (4.4d). We will use these values in the remaining analysis unless otherwise explained. 
We remind that we have mi^ < Mi6/smai. 



4.2 Lepton Flavour Violating Rare Decays 

Besides the neutrino mass measurements, the leptonic sector is already constrained from 
searches of LFV rare decay. In this specific model, the LEV couplings in Eqs. (2.5b) and 
(2.7) lead to LFV rare decays of charged leptons ia — )• i/3^- They are induced by similar 
loop diagrams as the one leading to neutrino masses. Using the general result for one loop 
LFV rare decays [16], we find for the decay width 
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where we have neglected the contributions from the suppressed U(l)x-breaking couplings. 
LFV rare decays are already strongly constrained by several measurements [12] to 

(4.11a) 
(4.11b) 
(4.11c) 

The MEG experiment aims at improving on the present bound down to Br(;U — )• 67) '^ 10"^'^ 
and has released the first result last year, Br(/i — )• 67) < 2.8 x 10^^^ at 90% C.L. [17]. The 
bounds on the LFV r decays will be further improved by a Super-B factory [18]. 

Inserting the values for the couplings in Eq. (4.7), we find that the bounds on Br(r — )• 67) 
and Br(r — )• ^7) can be readily satisfied even for values of itirr as small as 100 GeV and 
g^^r as large as 0.2. For ge,g^ ~ 0.1, the bound on Br(/i — ;• 67) points towards relatively 
large values of turr, turr > 1.1 TeV. However, taking 5^ ~ 0.02 and ge ~ 0.01, turr 
as small as 100 GeV can still be compatible with the present bound on /i — )• 67. Notice 
that, for such values, the expression for neutrino masses, Eq. (4.5a), implies g^ ~ O.Olg 
so despite relatively small g, the hierarchy imposed by the approximate U{l)x symmetry 
(i.e., ^A <C g) is still satisfied. An alternative possibility is ge <^ g^ or ge ^ g^. In the 
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case Qe "^ g^ii the ea elements of the neutrino mass matrix should be accounted for by the 
gg^r] contribution; i.e., {m^)ea = [ge{gA)a + ffa(5A)e]^- Similar consideration holds also 
for the case ge^ g^i- 

Since the processes /i — t- eee and fi — t- e conversion on nuclei cannot proceed at tree 
level in this model, the contributions in both cases are one-loop effects and are dominated 
by the effective vertex fla^'^ PieFfj^i, . The bound on this effective coupling from ^u — t- 67 is 
stronger than those from fj, — )• eee and /i — t- e conversion [19] because the former is a two 
body decay. 

4.3 Anomalous Magnetic Moment of the Muon 

Similarly to the LFV rare decay, the new couplings in Eqs. (2.5b) and (2.7) give rise to 
magnetic dipole moments of charged leptons £a 

where Xqq, is defined in Eq. (4.8). This contribution leads to a deviation of the dipole 
moment of the muon from the SM prediction 

fe, = Us - 2)„/2 . -< 1./ ~ 2.4 X 10- m^m , (4.13) 

where we have neglected the U(l)jic-breaking couplings. The present uncertainty on a^ is 
6 X 10"^'' [12] so the deviation is below the current experimental sensitivity and theoretical 
uncertainty. After an improvement on the theoretical and experimental uncertainties, the 
muon anomalous magnetic dipole moment, (5a^, will be a powerful test of this model. 



5. Dark Matter 

As discussed in sec. 2, the lightest neutral scalar, (5i, is stable due to the Z2 symmetry and 
is a candidate for the dark matter of the Universe. In order to be thermally produced in 
the Early Universe with the right abundance, the annihilation cross section needs to be 

(cr((5i(5i -^ anything)!)) ~ 3 x 10"^^ cm^/sec , (5.1) 

where v is the relative velocity. More precisely, as the mass splitting between the lightest 
and next-to-lightest scalar particles might be small, we have to take into account both 
particles during freeze-out. This requires a calculation of the self annihilation cross section 
of both b\ and b^ as well as the coannihilation cross section of b\-b2- b^ then decays into 
b\ , so the number of DM particles today is equal to the sum of the numbers of b\ and ^2 at 
the decoupling. In the following, we shall discuss these modes and determine the range of 
the parameters for which each mode will be relevant. We shall then discuss the possibility 
of direct and indirect detection of dark matter within the present model. 
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5.1 Dark Matter Abundance in the Universe 

Several processes contribute to the annihilation of dark matter in the Early Universe. They 
are depicted in Fig. 2. Their relative importance depends on the choice of parameters of the 
model. For the typical values of masses and the couplings that we consider, the dominant 
annihilation modes determining the dark matter abundance are mediated by the Higgs and 
the remaining modes are negligible. 

Higgs-mediated annihilation: An important annihilation channel is the annihila- 
tion via Higgs exchange. From Eqs. (2.3) and (2.5), we find that the coupling of 61^2 to the 
Higgs field is given by 

XiVHhdf = — ((Ahai - A//A2) sin^ ai + Xh,j, cos^ qi - 2Xha^ sin ai cos ai) h6f 

Ml - fil cos2 ai - n\ sin^ ai) ^^■'^^ 
^-h6f 



VH 

with i = 1, 2. For simplicity we neglect the subdominant U(l)x-violating terms. The 6i62h 
coupling is absent in the CP conserving Higgs potential. The CP violating terms lead to a 
coupling 6i52h and therefore induces coannihilations, which suppress the DM relic density. 
At low values of masses. Mi <^ rnw, with m\\r the VF-boson mass, the annihilation 
into fermion pair final states is important. The cross section for this channel [20] is given 

by 

{a{SiSi ^ ff)^v) = A^e^^ (^^,i^,^, ^^ , (5.3) 

where m/ is the fermion mass for the kinematically accessible channels and A'^c = 3 (1) for 
quarks (leptons). In the limit 5 <^ 2Mi, we have to take into account the annihilation of 
62S2 in the calculation of the DM abundance as 

(cr((52(52 -^ fDnv) ~ {a{5i5i -^ fDuv) ■ 

For Ml > mfj^r, the Higgs-mediated modes can dominate the freeze-out processes. For light 
dark matter, Mi < rrir, they can annihilate only into light fermions so the cross section 
is suppressed by the small fermion masses and cannot provide the dominant annihilation 
channel at freeze-out. 

For heavier masses. Mi > myy, three-body decays and decays into gauge bosons need 
to be taken into account. As shown in [21], even for 70 GeV < Mi < my/, the three body 
annihilation mode 5i5i — t- /i* — t- WW* — t- Wff is comparable to or can even dominate 
over 6161 —7- //. For Mi > mw, the Higgs-mediated DM annihilation into a W boson pair 
becomes kinematically allowed and soon dominates. Its annihilation cross section is given 

by 



{ai6i6i^WW)v)H- ^^_ 2 



9sU{2)\^l\ 

327rM2 




(5.4) 
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Figure 2: Different dark matter pair annihilation channels. 
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As the mass increases, more channels such as annihilation into Z boson pairs and top pairs 
become available but, the annihilation into a W boson pair will still dominate the cross 
section. Regardless of the final states, the cross section of any annihilation mode through 
s-channel Higgs exchange can be related to the corresponding Higgs decay rate by 

{ai6i6i ^h* ^ final state)/^f) = i2mhT{h -^ final state)) |^^^2Afi ~^]i^ UM'^^-m'^y ' 

'(5.5) 
Hence the importance of each channel can be obtained from the branching ratio of the 
corresponding Higgs decay channel by identifying the Higgs mass with the centre of mass 
energy -EcM = 2Mi. Let us take our typical value Mi = 70 GeV to evaluate A^. In this 
case, according to HDecay [22], the Higgs decay width is r|2Afi = 8.3 MeV. Taking into 
account the contribution of 82 for 6 <^ Mi, for nih = 120 GeV, the dark matter abundance 
implies a Higgs-DM coupling of Al ~ 0.07. We will use this value for reference in this 
study, unless otherwise indicated. As it can be seen from Eq. (5.5), larger Higgs masses 
ruh require a larger Higgs-DM coupling A^. 

Annihilation via gauge interactions: Another possibly relevant contribution is 
due to gauge interactions, which in principle involves two different processes. The process 
shown in Fig. 2(i) is suppressed by the heavy mass of A"*", but, for Mi > mw, the one 
shown in Fig. 2(h) results in the cross section 



-.4 / 2 




' (5.6) 

A comparison with the Higgs-mediated annihilation into W boson pairs shows that the 
annihilation via gauge interactions is subdominant for small mixing in the neutral sector 
as can be seen by computing the following ratio of cross sections 



{a{5i5i^WW)v)H _ \\lV ( vjj 



{a{5i5i -> WW)v) sin^ ai sin^ 02 v4Mf - m\ 



(5.7) 



This cross section can lead to the correct DM abundance for large mixing in the scalar 
sector, allowed for Mi > mz/2, and typical masses in the range Mi ~ 100 GeV — 200 GeV. 
We do not consider this possibility further. The three body final state annihilation pro- 
cesses, see Fig. 2(k), are related in a similar way. If the annihilation via gauge interaction 
dominates, the box diagram with two W bosons can be significant in the DM nucleon 
interaction and therefore in the direct DM detection experiments. 

Annihilation into Higgs pair: For Mi > rrih, the annihilation into a Higgs pair can 
take place and even dominate over other channels. In particular, the quartic couplings of 
H, $ and A can lead to 
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For simplicity, in the following we consider only the region Mi < rrih, where this annihila- 
tion channel is absent. 

Annihilation into neutrino and anti-neutrino pairs: A dark matter pair can 
annihilate into pairs of (anti-)neutrinos through the t-channel chirality-flipping diagram 
shown in Fig. 2(b). To leading order in U(l)x-violating couplings, the cross section is 
given by 

{a{6i6j -^ Uai^fj)v) = {a{5i6j -^ h'aVp)v) = , . ^ |5a(5A)/3 + 9p{9d.)a? (5.9) 

with i, j = 1, 2 in the limit Mf ^ ^rr- 

For typical values, the cross section can be estimated to be 

,,,, ,, ^ ^^ 07 /300GeV\ Vsin2 2ai\ / ^A \V 9 \2cm3 

< {a{5i5i -^ anything)^) (5.10) 

and it cannot therefore explain the DM abundance. 

A quantitative connection between neutrino masses and the dark matter abundance [2] 
can be obtained if a different appropriate choice for the U(l) symmetry, protecting neutrino 
masses, is made. In this case the dominant dark matter annihilation channel can be 
6161 — )• vu, vv. This leads to an upper bound of the order of 300 GeV on the masses of -E]^, 
^'r 1 ^R c^iid z/^ making their production at the LHC possible. 

Annihilation into lepton pairs via heavy fermion exchange: In addition, pro- 
cesses such as 6i6j — t- l^l^,ui> can also take place via heavy- fermion exchange, which is 
shown in Fig. 2(c). The cross section can be written as 

to leading order in the final state lepton masses. The p-wave contribution vanishes. The 
estimate for the dominant annihilation into r leptons with our typical values results in 



Qt 



0.1 



4 (702 + 3002)2GeV^ cm' 



a{5,5, ^ T-T+)v) ^ 3.3 x 10^^^ -^ \,., ,^ , 2 n2 TIT ' (5-12) 



{MiMj + m\j^Y sec 



which cannot give a dominant contribution to the dark matter abundance. Note that this 
cross section depends on some of the parameters which control neutrino masses and the 
rate of LFV processes and cannot therefore be enhanced arbitrarily by a different choice 
of parameters. 

DM annihilation into photons: There are several one-loop diagrams that con- 
tribute to annihilation into a photon pair, as summarised in Fig. 2(e). For the values of 
Ml up to the electroweak scale, the cross section can be estimated as 

,/,-,- s > a^Gl.sin'^Qi r, n. /sinai\ / Mi \ cm^ 

W<SA^77V> 5if^*^N3xlO-«(^) (^^) — , (5.13) 
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which is neghgible. 

Coannihilations of S\ and ^2 via the Z boson: Through the mixing with the 
neutral components of the triplet A, 6i and 82 couple also to the Z boson (see Eq. (3.6)). 
This coupling allows the 6162 coannihilation into kinematically allowed modes, such as vu, 
e~e~^ depending on the values of Mi and M2, shown in Fig. 2(d). The cross section can be 
evaluated as 



/ ^x X . ff\ \ .. G^, sin^ ai sin^ 02 
[a{6i62 -^ ff)v) = Nc—^ 



27r 






m^y 



(5.14) 



where a\ + a\ is given by 

-cos^ 29w + s\ii^ Ow ■, T, — (9-4cos26'vK + 4cos46'h/) , 

— (6 + 2 cos 261^/ + cos 461,4/) , (5.15) 
36 

for leptons, neutrinos, up-type quarks and down-type quarks. Hence a\ + a?^^ 0.13 — 0.25. 
If the mass of the dark matter is larger than that of the h quark but smaller than mz, 
this mode of annihilation can be estimated for our typical values in Eq. (3.9) to be 

/ / ,- ,- -X \ / sin ai sin tto \ / qv 6"^ no 9\ cm^ 

(aiiA - bb).> =. [-^j^) (7 X 10-- j^^^^^j^ + 4 X 10-- v^) — . (5.16) 

where the p-wave contribution dominates over the s-wave for typical value of u ~ y^1/20 
at the DM decoupling time. For our typical values, this is only about < 1% of the total 
annihilation channel cross section (Jj(5j — )• anything). Unless |sinaj| > 0.5, we expect 
coannihilations not to significantly modify the amount of dark matter at freeze-out with 
respect to the case in which they are neglected. Notice that in the case of Mi ~ 7tt-^/2, there 
is an enhancement of the coannihilation cross section, which can contribute significantly 
to freeze out. 

^2 decay: 82 eventually decays into 5i and neutrinos or other kinematically allowed 
light fermions via Z exchange. It can also decay into lepton-pairs via the R — R' exchange. 
To leading order in the U(l)x-breaking couplings and the mass splitting (5, the decay rates 
are 

T-/X R. X \ rfx ^. X --\ 5^sm^2ai _ 2 /c1v^ 

T{d2 > diuv) = T[62 > divv) ~ n 3..r2 2 baffA/j + fi'/sSAal , (5.17a) 

1920(1 + dai3)T^ Mfmjjj^ 

r{62 A 6iuu) ~ ^^ sin^ ai . (5.17b) 

For 6 S> iTie, the decay rate into electron final state pairs is 

T{62 A <5ie+e-) ~ 4(ai + a^) T{62 A Siui?) = (cos^ 2ew + 4 sin^ 9w) ^{^2 ^ ^J^i^) 
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This channel is suppressed with respect to the corresponding neutrino channel by a factor 
of 2. The same decay mediated via the new heavy fermions has a rate given by 

Ti6, ^5,e e , 6,.u) :. ^,,^s ^^2 _ ^2^^y ■ (5-19) 

For the typical values we have chosen for the scalar and leptonic parameters, the decay via 
the Z boson into a neutrino anti-neutrino pair dominates with a decay width 

Hence, the decay happens before big bang nucleosynthesis (BBN) and does not affect the 
BBN predictions. We comment on the case with a small mass splitting 6 ~ 0(10 — 100) keV, 
as it is required in the low mass region to explain DAMA via inelastic scattering, in the 
next section. 

5.2 Direct Dark Matter Searches 

Direct DM searches look for interactions of DM with the nuclei (electrons) in the detector. 
The differential scattering rate is given by (see e.g. [23]) 

^(En,t) = -^,[U/fnZ + {A-Z)fa^F\E^) f''^ ^^JjSS^!^^M. , (5.21) 

where Er is the recoil energy of the nucleus, p^ is the local DM density. Mi is the DM 
mass, nir is the reduced mass of the dark matter-nucleus system; v is the speed of dark 
matter relative to the nucleus; fp/ fn is the ratio of the coupling of DM to protons compared 
to neutrons, F{Eji) is a nuclear form factor describing the nuclear structure, fiocai is the 
local DM velocity distribution, Vesc is the escape velocity and finally (T„ is the dark matter- 
neutron cross section. 

an and fp/fn depend on the dominant DM-nucleon interaction which can be obtained 
by rotating the diagrams in Fig. 2(a) and Fig. 2(d). There is also a contribution from two 
W boson exchange box diagram, which is negligible for our typical values, but becomes 
important above the W boson mass threshold and large values of sinai. The t-channel 
Higgs boson exchange leads to [20] 

(5.22) 
with fp/fn ~ 1- Here p,Sin is the reduced mass of the dark matter-neutron system, irip is the 
nucleon mass and / parametrises the nuclear matrix element, 0.14 < / < 0.66 in [20]. Note 
that the Higgs-mediated cross section strongly depends on the uncertainties in the nuclear 
matrix element. This interaction would lead to elastic spin-independent (eSI) scattering 
and to nuclear recoils which have been extensively searched for by various experiments. 

If the mass splitting between 61 and 62 is small (0(10 — 100) keV), Jin — ;• 62n through 
t-channel Z boson exchange becomes kinematically allowed. In the limit 5 <^ fJ-Smi the 



17 



cross section is given by 

sin^ ai sin^ a2gFj»lin ?^ 1.3 x 10"^^ / sinaisina2 \ ^^2^ ^^^S) 



Cn = 



with fp/fn = —(1 — 4: sin 9\y) ps —0.08. The dependence on 5 is reflected in the lower 
limit of the integral in Eq. (5.21). This accommodates the inelastic spin-independent (iSI) 
scattering scenario, in which dark matter, 5i, is converted into a slightly heavier particle, 
62, while scattering off nuclei [24] (see also [25, 26]). It should be pointed out that our 
model accommodates the small mass splittings required, 6 ~ (m^ + sin a im^^) /Mi, as it 
is naturally suppressed by the U(l)x-breaking terms. 

Experimental constraints:: In the following, we discuss the experimental con- 
straints from direct searches ^. Many experiments have searched for nuclear recoil signals, 
e.g. XENONIO [29], ZEPLIN-III [30], CRESST-II [31], KIMS [32] as well as PICASSO [33], 
and, recently, the CDMS II [6], the CoGeNT [8] and the XENONIOO [9] experiments. The 
most stringent bounds on the spin-independent elastic cross section come from CDMS-II 
and XENONIOO: for Mi = 55 GeV o-„ < 3.4 x IQ-^'^cm^ at 90% C.L. from the XENONIOO 
first data and for Mi = 70 GeV Un < 3.8 x lO'^^cm^ from CDMS-II. Evidence and hints 
of dark matter detection have also been reported but await further confirmation. The 
DAMA/LIBRA experiment in Gran Sasso searches for an annual modulation of the DM 
scattering signal due to the Earth orbit around the Sun and the consequently annual change 
in the DM velocity relative to the detector. It has been reporting a positive signal for 13 
years [7]. The effect which is seen by DAMA at 8.2 a is refuted by other experiments 
attempting to directly detect the dark matter. The CDMS-II experiment reported two 
candidate events requiring a la-allowed region in the Mi — an plane roughly between 
21 GeV < Ml < 51 GeV and an — 10~^^ cm^-10~^^ cm^ for eSI scattering being consistent 
with all other null results. The allowed values of an extend up to (7„ ~ 10^"^^ cm^ for low 
masses if the bounds from XENONIOO can be relaxed [9]. The CoGeNT experiment sees 
an excess of events at very low energies below 3 keV, which, if not due to backgrounds, 
can be interpreted as dark matter- nucleon eSI scattering with Mi ~ 7 GeV-11 GeV and 
ajy ~ 3 X 10~^^ cm^-1 x 10"^'' cm^ (for other analysis of this and other DM direct searches 
data, see also Refs. [34, 35]). These results can be compatible with the DAMA preferred 
region for an intermediate amount of channelling but are in tension with the XENON and 
CDMS results. Recent analyses of the relevant experiments have been performed [34] (see 
also [35, 36]) in order to obtain global limits on the DM- nucleon cross section and indicate a 
tension between DAMA, CDMS, XENONIOO and CoGeNT data. However, the results still 
depend strongly on the underlying assumptions of the experiments, which are not settled 
yet. For example, a different choice for the effective light yield of the XENONIOO [9, 11] or 
the channelling in the DAMA experiment can lead to significantly different allowed region 
of parameter space. Taking a conservative effective light yield for the XENON experiment. 



^The proposed explanation of the DAMA signal by scattering off atomic electrons [27] is disfavoured by 
different analysis due to the tension between the DAMA spectral data and the modulated signal [26, 28] as 
well as the loop induced interactions with nuclei [28]. 
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the combination of Fig. 2 in [34] and Fig. 3 in [20] (also in Fig. 2 in [35] and Fig. 5 of 
[36]) suggests that there might remain a region of parameter space ~ lOGeV explaining 
DAMA or CoGeNT by eSI scattering via Higgs exchange which is compatible with the 
bounds from other experiments. However, this region is excluded if the effective light yield 
of XENONIOO is higher or the region allowed by DAMA is more restricted. The cross sec- 
tion required by DAMA and CoGeNT needs an intermediate amount of channelling and/or 
a sizable source of background at low energy in CoGeNT. A stronger tension between the 
regions preferred by DAMA and CoGeNT with the one required to explain the two CDMS 
events remains and could be partially alleviated only by assuming a different DM velocity 
distribution (see e.g. [35]). We do not therefore restrict ourselves to one analysis but base 
our discussion on the analyses in [37, 38] and [34, 35, 36]. 

Elastic DM-nucleon scattering: The DM-nucleon cross-section in Eq. (5.22) is 

controlled by the same parameters as the dominant annihilation cross section, Eq. (5.5). 
After factoring out Al, we find 



/Vl n^p(4^i - "^D^ HSi6i -^h* ^ SM final states)^; ) 
irMimj^vfj 4T{h — )• SM final states) |m^_5.2Mi 



a„(~ dp) = _„^_4.,2 .T^^i. , cA^fi,.„i^.^.„^M :— ■ (5-24) 



Thus, once the DM mass is fixed, the DM-nucleon cross section is uniquely determined and 
depends mildly on the value of the Higgs mass. Ignoring the positive signal in favour of DM 
scattering, the recent analysis leads to a typical cross section o"„ < 2 — 3 x 10~^^cm^ for a 
DM mass of ~ 10 — 130 GeV. These values can be accommodated in our model, depending 
on the value of the Higgs, DM mass and the nuclear matrix element. We expect soon 
a positive signal in direct detection experiments unless Mi ~ m/i/2 (see Eq. (5.24)). At 
Ml ~ 70 GeV, CDMS and XENONIOO are already constraining the parameter space and 
large values of /, / > 0.2, are not compatible with the bounds from direct DM detection 
experiments. For smaller values of the mass, the elastic cross section is suppressed by 
the cancellation between the Higgs mass and 2Mi, see Eq. (5.24). For example, for a 
slightly smaller DM mass, e.g. Mi = 65 GeV, the corresponding total Higgs decay width 
r = 5 MeV leads to Xl ~ 0.04 and therefore a value of the cross section a"„ ~ 1.8 x 10~^^cm^ 
for / = 0.3 which is below present bounds. 

If DM has been observed and the two candidate events of CDMS-II are due to dark 
matter, the allowed region in the parameter space. Mi '^ 20-50 GeV and o"„ ~ 10^'^^ cm^- 
10~^^ cm^, can be accommodated within our model via scattering by Higgs exchange with 
e.g. Ml = 50 GeV and a light Higgs ruh = 120 GeV, and (T„ ~ 5.4 x lO"'^^ cm^, where 
we have used the lower bound on / = 0.14. In the case of sizable coannihilations, the 
required DM coannihilation cross section is smaller and therefore the elastic DM - nucleon 
cross section is reduced. This improves the consistency with the two candidate events of 
CDMS-II. The recent CoGeNT results, if interpreted as a dark matter signal, require a 
different region in the parameter space, with smaller masses 7 GeV < Mi < 11 GeV and 
higher cross sections (T„ ~ 10~^^ cm^-lO"^" cm^, which might also explain DAMA for 
intermediate channelling, as discussed above. In our model, for fixed Mi in the range of 
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interest we predict the value of the cross section 



in agreement with the experimental results for small Mi and larger /. It is curious to 
notice that this is exactly the range where eSI solution with channelling for DAMA comes 
close to the preliminary XENONIOO bounds [9]. On the other hand as shown in [10], in 
this range the bounds are sensitive to astrophysical uncertainties such as the dark matter 
escape velocity. In future, more robust bounds might conclusively refute this solution for 
DAMA. In this case, our model is still compatible, as Mi can take on higher values. 

Inelastic dark matter: The inelastic SI scenario has recently attracted much interest 
as it can simultaneously accommodate the DAMA signal and the CDMS data [6] because 
scattering off heavy nuclei (such as ^■^ I) is favoured with respect to the one onto light 
nuclei and the modulated signal is enhanced compared to the unmodulated one. The 
latest global analysis [34] finds three possible regions for Mi = 10 GeV, 40 GeV, 50 GeV. 
However, the regions around Mi = 40 GeV, 50 GeV are both excluded by the CRESST-II 
data as discussed in [34] and by the bound of Super-Kamiokande on the neutrino flux from 
DM annihilations in the Sun [39] . We will focus on the lowest allowed values of the cross 
section in the region around Mi = 10 GeV: dp ~ 1 x 10"^" cm^, which roughly corresponds 
to an = 3.3 X 10~^'^cm^. This region is also suggested by [38]. The DM-neutron cross 
section can be estimated from Eq. (5.23) as cr„ ss 6 x 10~^° (sinai sin 02/0. 07) cm^, where 
we have used the largest allowed value of sinai, close to the present bound from the 
invisible Z-decay and testable with a moderate improvement of these searches. 

A comparison with Fig. 7 of [38] shows that the resulting cross section can still explain 
DAMA, whereas the analysis [34] already excludes this cross section assuming a standard 
DM profile. In general, iSI scattering strongly depends on the velocity distribution because 
only the high energy tail can scatter. Hence, this region is probably still allowed due to 
the astrophysical uncertainties. Concerning the values of masses, it requires fine tuning at 
the % level to obtain a light DM mass Mi in the region of interest. As already mentioned, 
the mass splitting (5 ~ 20 keV is naturally small due to the V{l)x symmetry and could be 
obtained for example for m^,7fi^^ sinai ~ 2 x 10~^ GeV^. In this case neutrino masses 
would be dominated by the contribution due to ggi^rj, as rji'^^^^AA) i^ *oo small. For 
a small mass splitting, 5 ~ 0(10 — 100) keV, as it is required in the low mass region to 
explain DAMA, the decay of 62, which survives after freeze-out, happens at very late 
times and after galaxies have formed, as r^j ~ 7 x 10^^ (20 keV/5) (0.07/sinQi sin 02) s. 
Its effect can be estimated by looking at the energy densities pi = pi + P2 — Pi, Pi denoting 
the energy density of 6i. Since 6i are non-relativistic, the energy densities are given by 
Pi = riiMi. The decay leads to a negligible energy density increase for neutrinos 

PV PV iV 2Mi \ly 

where ^y > Y^A?7i^j,^/91.5eV rs 5 x 10~^ has been used as lower bound for fij/. As their 
energy is too small E^ < 5 ^ 20keV, they evade detection in neutrino detectors and. 
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since the energy of neutrinos is below the nuclear binding energy, they cannot destroy the 
outcome of big bang nucleosynthesis. 

Notice that, in addition to the iSI scattering, elastic scattering will necessarily be 
induced with a cross section determined by Mi as discussed above. We expect a large cross 
section for elastic scattering in addition to the inelastic one, which can be compatible with 
present bounds from CDMS and XENONIO for small values of / and/or a conservative 
treatment of experimental uncertainties [40] and astrophysical parameters [10]. 

If the mass splitting 6 is even smaller, the exothermic dark matter (exoDM) [41] 
scenario might explain DAMA and the other direct detection experiments within our model. 

In summary, present direct dark matter experiments provide contradictory results, 
with DAMA showing a strong evidence of annual modulation of the signal, CDMS and 
CoGeNT showing possible hints in favour of DM if their signal is not due to backgrounds, 
and the other experiments reporting null results in the region of the parameter space 
of interest. If we dismiss the possible positive signals found so far, our model typically 
predicts an elastic cross section within the reach of present and future experiments, unless 
Ml — )■ ?7i/i/2. Otherwise, if we take the positive signals as a direct observation of dark 
matter, various possible explanation can be accommodated in our model, depending on the 
values of the parameters. For Mi ~ (21 — 51) GeV we can explain the CDMS two-events, 
while DAMA, with intermediate channelling, and CoGeNT require much smaller masses. 
Ml ~ (7 — 11) GeV, and a correspondingly higher cross section. The latter signals can 
also be explained with the iSI scattering which requires a value of sin ai close to the upper 
bound, sin^ ai ~ 0.07, and therefore testable in the future in invisible Z-decay searches. 
In the near future, new results for direct DM search are expected and in particular further 
data from XENONIOO experiment [42] and CRESST will help to clarify these issues. 

5.3 Indirect Dark Matter Searches 

Indirect dark matter searches look for gamma-rays, neutrinos, positrons, anti-protons and 
anti-deuterons from the regions of the galaxy or astrophysical objects (in the case of neu- 
trinos) where the concentration of dark matter is expected to be relatively high and anni- 
hilations are therefore strongly enhanced. A study in Ref. [43] has recently derived limits 
on the different dark matter annihilation channels leading to electron positron production 
by studying radio and gamma ray from galactic center. Moreover, in Ref. [44] bounds 
on different DM annihilation modes have been derived from Fermi-LAT diffuse gamma 
ray data. A comparison of the different annihilation channels with Fig. 2 of [44] shows 
that our model is not constrained by the Fermi-LAT data but a future improvement on 
the sensitivity will be able to provide useful constraints for light ©(few GeV) dark matter 
masses, when the dominant annihilation is into light quarks or rs. The coannihilation of 
DM into photons at one loop which has been estimated in Eq. (5.13) is well below the cur- 
rent bounds from EGRET and Fermi-LAT [45], too. The anti-deuteron cosmic ray search 
experiments AMS-02 and GAPS can test the DM annihilation to hadronic final states in 
the region around 0(100 GeV) [46]. 

Dark matter can also be captured in compact objects such as the Earth and the Sun 
due to the scattering on nuclei. This leads to a large flux of neutrinos either prompt 
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from the annihilations or as subsequent decay products from annihilations into charged 
leptons and quarks. This feature can be tested in present and future neutrino detectors 
such as SuperKamiokande and IceCube. These detectors will measure the total neutrino 
flux and can in principle determine the neutrino spectrum, if a sufficient energy resolution 
is available [47]. Notice that, although the overall detection threshold of IceCube is rather 
high, its DeepCore component has a threshold of 10 GeV [48] and can be used for this 
purpose, if Mi > 10 GeV. If the nucleon-dark matter interaction, as well as the dark 
matter annihilation, dominantly proceed via Higgs-exchange, the effect can be significant. 
In this case, for masses below 70 GeV, the dominant annihilation channel is into 6-quarks as 
the cross section scales with the final fermion mass squared. We therefore expect a rather 
soft neutrino spectrum with a fixed branching ratio into c-quarks and rs, below present 
constraints [49]. For heavier masses new channels are open: annihilations into gauge and 
Higgs bosons lead to a hard neutrino spectrum which can be more easily detected at present 
and future detectors [50, 51]. In the iSI case, with a mass splitting in the 0(10 — 100) keV 
region, due to the high inelastic scattering cross section, dark matter would be copiously 
captured in the Sun [39, 52, 51]. A population of 82 particles would form, which could 
subsequently decay into 5i along with low energy neutrinos, not detectable with present 
techniques, or could annihilate as discussed above. For the inelastic scattering cross sections 
and required DM mass Mi (see the previous section), the hard channels, as annihilations 
into rs and into neutrinos, are already constrained by Super-Kamiokande data to give a 
subdominant contribution, but annihilations into bs and cs are allowed [39] . We recall that 
in our model, the main annihilation modes are ^i^i — t- /i* — t- 66 and ^2^2 — s- /i* — )■ 66 giving 
a rather soft neutrino spectrum. Therefore, it is possible to explain DAMA with inelastic 
dark matter evading the present constraints from dark matter neutrino searches from the 
Sun. 

6. Other Constraints on the Model and Laboratory Signatures 

In this section we discuss electroweak precision observables which might constrain the 
model further and speculate about possible collider signatures. 

6.1 Electroweak Precision Tests 

As the additional particles are close to the electroweak scale and are charged under the SM 
gauge group, they lead to corrections to the electroweak precision parameters [53, 54, 55]. 
It has been pointed out in [55] that all contributions of physics coupling only to the lepton 
sector can be condensed into seven effective oblique parameters. The dominant effects are 
contained in the quantities S*, T, W^ Y [54, 55]. 

A study of an additional vector-like lepton doublet [56] shows that the contributions 
to S and T parameter exactly cancel out, because R and R' have equal masses, while W 
and Y receive tiny corrections 

^^Shmrr^ and F = "^4^ > 
1207r2 m\^ 1207r2 m\^ 
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respectively. 

We can neglect the contribution of (j) to the electroweak precision parameters because 
it is suppressed by a factor of | sinai sina2| relative to that of A. For the latter, the direct 
calculation of the wave function renormalisation results in 



247r2^' 5767r2 m^^^ ' 7207r2 m\' 4807r2 m^ 

where the relation 2 m^+ = m^ + m^_,.+ has been used and the results have been expanded 



m 

?^^=^iV^ = A„.,4. (6.2) 

It can be easily seen that the two additional fermionic doublets with opposite hypercharge 
as well as the triplet without VEV have a well defined decoupling limit. ^ can be chosen 
such that it cancels the contribution from the SM Higgs, relaxing the upper bound from 
the electroweak precision data on the Higgs mass. Without cancellation (i.e., for a light 
Higgs mass), the T parameter constrains ^ ^ 0.1 which translates into a bound on the 
splitting of the components of the triplet. This results in a mild bound on XhA2, e-5-i for 
m/\ ~ 500 GeV, the bound is Xha2 ^0.5. The other electroweak precision constraints are 
readily satisfied. 

6.2 Signatures at Colliders 

Higgs Boson Searches: If Mi < rah/2, the coupling which is responsible for the DM anni- 
hilation also leads to the decay of the SM Higgs boson into DM particles. For Az, > mi,/vH, 
its branching ratio becomes significant and even dominates over the decay into bb. This hap- 
pens for our typical parameter set, where we have A^ ~ 0.07. Hence, a light {mh < 2m\Y) 
SM Higgs decays dominantly into 5i5i or 6262- The DM particles 61 escape the detector. In 
case that the mass splitting between 62 and 61 is less than twice the electron mass, 62 will 
decay only into 5i and neutrinos which are also invisible. For larger mass splittings, the 62 
decay into e~e~^ can lead to a displaced vertex, which opens a new and distinct channel for 
discovering the Higgs [57], provided that the decay takes place inside the detector. For this 
to happen, the diameter of the detector, d, bounds the decay width F^j by dTs^/^j > v 
with V being the velocity of the particle 62 and 7 = (1 — t;2)~^/2. Assuming a dominant 
decay via the Z boson, this translates into a bound on the mass splitting 5 

5^in^ai>607r3^. (6.3) 

Hence, for the maximally allowed mixing sin2 ai ~ 0.07 in the case of Mi + M2 < mz, the 
ATLAS Muon detector [58] with a diameter of 22 m already requires a mass splitting of 
6 > 480MeV(7u)^'^. This displaced vertex would be a clear signal for a neutral next-to- 
lightest particle with SM couplings. 

Prospects for the LHC: Since this model contains several particles with masses in 
the reach of the LHC, we expect a rich phenomenology within the upcoming years. If the 
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new particles are not too heavy, the charged particles A++, A+, E^ and E'^ as well as 
the neutral particles 83, 5^, ur and u'^ can be produced through electroweak interactions. 
They will then decay into the SM particles plus 5i or 82- At the LHC, 5i appears as a 
missing energy signal. 82 subsequently decays into 8ivi> or, if kinematically possible, into 
8ie~e'^ . If the decay happens outside the detector or is into neutrinos, the displaced vertex 
cannot be observed and this decay will contribute to the missing energy signal. Since the 
masses of the components of the electroweak triplet A fulfil the relation 

2m\+ = m\++ +m\ , 

it might be discovered by measuring the masses of A+, A^^ and 8^, at the LHC, as long as 
^HA<f) and XhA(j> are small ttt-a — M3 ~ M4 (see Eq. (A. 4)). In fact, from the electroweak 
precision data. A"*" and A"'"'' are expected to be quasi-degenerate with a small mass splitting 
of m?^++ — "1^+ = XhA2vJj/'2- The coupling Qa can be determined by measuring the decay 
modes of E~^ because the branching ratio Br(£'^ — )• ^~<5i,2) oc \ga\'^ ■ The values of the 
components of g^ can be derived from a study of the decay modes of A"*" and A'''^. In 
particular, r(A++ -^ C^^'^1,2) « \{9A)agi3 + (^AJ/SffaP- By directly extracting g and (Ja 
at the LHC, it will be possible to cross-check the information on them from rare decays 
and the neutrino mass matrix (see sec. 4). 

As discussed in sec. 3, an improvement of the uncertainty on the invisible decay width 
of Z can test the model for Mi < mz/'^- LHC, being a Z factory, can in principle improve 
the precision of the r(Z — )■ invisible) measurement. 

7. Conclusions 

In this paper, we have proposed a model that simultaneously explains the missing mass 
problem of the universe and the tiny neutrino masses. In addition to the SM particle con- 
tent, there are only a complex scalar singlet and triplet as well as a vector- like electroweak 
fermionic doublet. We impose an approximate \]{l)x symmetry which is broken to a rem- 
nant Z2 symmetry. The unbroken Z2 symmetry guarantees the stability of the lightest 
scalar in the model, 81, a quasi-singlet of SU(2)x,, which plays the role of dark matter. 
In the limit of exact U(l)x symmetry, neutrinos are massless and only after the U(l)x is 
broken to the Z2 symmetry, neutrinos acquire a mass term at the one-loop level. The Z2 
symmetry forbids a tree-level neutrino mass term and the usual seesaw mechanism does 
not take place. Hence, the smallness of neutrino masses is explained by the small breaking 
of the U(l)x symmetry as well as the loop suppression. With the minimal particle content 
of the model, one of the neutrino mass eigenvalues vanishes and the neutrino mass scheme 
is therefore hierarchical. In order to obtain a non-hierarchical neutrino mass scheme, the 
minimality of the model has to be relaxed and more vector-like fermionic doublets have to 
be added. The strongest constraints come from searches for lepton flavour violating pro- 
cesses, in particular fi — )• ej, which already probes the relevant parameter space. Future 
searches for /i — t- 67 will provide a very sensitive test of our model. 

In this model, DM is produced thermally in the Early Universe. We discussed the 
different dark matter annihilation channels and identified the dominant one to be the one 
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via Higgs exchange, for Mi <C rriw. All other channels are sub dominant. The predicted 
cross section is compatible with the value required to explain the observed DM abundance. 

The interactions responsible for DM freeze-out induce also scattering of dark matter 
off nuclei, relevant for direct DM searches. For our typical values Mi = 70 GeV, 6i scatters 
elastically via Higgs exchange. The obtained scattering cross section is just below the 
current experimental bound and moderate improvements on the sensitivity can probe part 
of the relevant parameter space. Our model can also accommodate light dark matter 
with mass in the few GeV range, which has been invoked to explain the CoGeNT and 
DAMA results via elastic scattering. This process is mediated by the Higgs exchange 
and can have the required value for the cross section. For heavier masses, Mi ~ 20- 
50 GeV, the two events recently reported by CDMS can be interpreted as dark matter 
elastic scattering with a cross section which is compatible with the predictions of our 
model. The first results of the XENONIOO experiment [9] disfavour most of the parameter 
region of DAMA, CoGeNT and the two events from CDMS depending on the assumptions 
on astrophysical uncertainties [10] and the ratio between electron equivalent energy and 
nuclear recoil energy Cqq [11]. Further data from the XENONIOO experiment as well as 
other experiments is needed to resolve this uncertainty. We also studied the possibility of 
inelastic spin independent solution for DAMA. For small mass splittings and small dark 
matter masses, 6i can scatter inelastically to 62 via Z boson exchange through mixing 
between scalar singlet and triplet. In our model the mass splitting 6 can be naturally small 
due to the U(l)x symmetry. In order to accommodate the solution, the singlet-triplet 
mixing has to be relatively large and just below the upper bound from the invisible Z 
boson decay width. Thus, a slight improvement on the precision of the invisible Z decay 
width can probe this phenomenologically interesting part of the parameter space in our 
model. 

We demonstrated that bounds on electroweak precision observables do not constrain 
the model further. Even more, the upper bound on the SM Higgs from electroweak precision 
data can be relaxed. The new particles can in principle be produced at the LHC and will 
eventually decay into stable 61 which escapes detection. The second lightest scalar, 62, 
dominantly decays via Z exchange and might lead to a displaced vertex in the detector 
for sufficiently large mass splitting 5, or can decay outside the detector contributing to 
the missing energy signal. It is possible that H — )• 61S1 and H — )• ^2^2 dominate over 
the SM mode H — )■ bb, if ^1^2 are sufficiently light. In this case, the Higgs would decay 
mainly invisibly. The relevant coupling A^ for this decay is fixed by the DM annihilation 
rate. Collider searches for the fermionic doublet can be also performed. By studying the 
subsequent decay of the charged components of the doublet into a charged lepton, we 
can determine the Yukawa couplings of these particles to different fiavours. The fiavour 
structure of these couplings also determines the flavour structure of the neutrino mass 
matrix so this provides another method to cross check the model. 

Let us finally, comment on an alternative possibility, which leads to a tight connection 
between neutrino masses and the dark matter abundance [2]. If the guiding symmetry 
is not U(l)j(:, but an approximate lepton number U(l)i or \J{1)b-l, the dominant dark 
matter annihilation channel may be 5i5i — t- i'a'^/3,i^ai^i3 resulting in a direct connection 
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between the dark matter abundance and neutrino masses. This leads to an upper bound of 
the order of 300 GeV on the masses of E^ , E'^ , vr and v'j^ guaranteeing their production at 
the LHC. In this case, the neutrino flux from dark matter annihilations inside the Sun will 
be monochromatic with a general flavour composition determined by the flavour structure 
of the new Yukawa couplings of the model. As recently shown in [59], this can lead to a 
novel seasonal variation in IceCube which cannot take place in models predicting only a 
continuous spectrum or democratic neutrino flavour composition. 

In summary, we have presented here a model which explains simultaneously the origin 
of neutrino masses and the dark matter. A global U(l)x symmetry, explicitly broken to a 
residual Z2 guarantees the smallness of neutrino masses, generated at the loop- level, and 
the stability of dark matter. The model has a very rich phenomenology, such as lepton 
flavour violating processes, invisible decays of the Z-boson, collider signatures, which will 
make the model testable in the near future. Dark matter annihilations dominantly proceed 
via Higgs-exchange. Elastic and/or inelastic scattering off-nuclei can also be induced by 
the Higgs or Z exchange and can explain the possible signal or hints for dark matter 
direct detection which have been recently reported. So far, we have considered an explicit 
breaking of the additional U(l)x symmetry, but a version with a gauged U(l)x symmetry 
is in preparation. 
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A. Scalar Mass Spectrum 

The terms in Eqs. (2.3) and (2.6) with the vacuum expectation values defined in Eq. (3.1) 
lead to the following charged scalar masses 



2 



■^ 2 , ^HAl + A/f A2 2 ^ A 1 ^ 

m^++ =fi^ + Vh , (A. la) 

ml, =/.i + ^vj, . (A.lb) 

In order to obtain the mass eigenvalues of the neutral scalars, one has to diagonalise their 
mass matrix. Remember that we have decomposed A^ and (p as A*^ = (Ai + iA2)/\/2 and 
(f) = {(pi + i(/)2)/\/2. In the basis {(pi, (p2, Ai, A2), the mass matrix is given by 



m? 



/"^0i ml^+rhl^ \ 

ml 

ml ) 
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(A.2) 



where 



4 - 2A' + (Ah<^ - 2\h4) '^ = ml + ml, (A.3b) 



J2. 

m-^i = /^0 + 2/i^ + [\h<i> + 2Ah</,J ^ = m^ - m^ , (A.3a) 

^2 _ ,,2 or,2 , /\.. . o\.. A "h _ _2 , ^2 

"^A = /"A + (AhAI - A//A2) ^ , (A.3c) 

„2 

"^^A = -Aha^y" ' (A. 3d) 

^2 

"^^A = -A^Afli^ ■ (A.3e) 

The diagonahsation by a transformation into the mass basis given in Eq. (3.5) yields the 
mass eigenvalues 

(A.4a) 



Ml = \ (ml, + ml - ^(ml - ml,? + 4 (m^^ - rf.^^)^ :. r.^, - ^"^^^ _^f ^^ 



(A.4b) 



,2 I ~2 \2 



Mi = I l^ml, + ml + ^{ml - ml^Y + 4 (m^^ + ml^f] ^ ml + ^^^^^f^^ 



(A.4c) 



Ml = J ( m^, +ml + \l{ml-ml,Y + ^(m--mlX\ :. ml + ^"''^ ^^^^' 



- 1 ,10^2 -r- '"A -r Y V"'A " "V2-' ^ "* ^^"VA ~ "VA^ I "^ '"A t ^2 _ 2 

\ / A <7i2 



(A.4d) 
where in the last equation, we have assumed 

"^A > rnl^^ml^ and m^ - m\^,ml - ml, > m^^ ± m^^ . 
Positiveness of M^^ guarantees ((/>) = (A) = 0. The mixing angles are 

smaicostti = = , (A. 5a 

2 / ^2 



[ml-ml^^ +4(m2^ + m2 



()iA 
sma2COsa2 = — , ^ = ■ (A.5bj 



+ 4 mL -mi 



'^A~'^<t>2j "T-* \^"VA ~ "VA 
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